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We consider the use of N spin-i particles for indicating a 
direction in space. If N > 2, their optimal state is entangled. 
For large N , the mean square error decreases as iV~ 2 (rather 
than N" 1 for parallel spins). 
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Information theory usually deals with the transmis- 
sion of a sequence of discrete symbols, such as and 1. 
Even if the information to be transmitted is of continu- 
ous nature, such as the position of a particle, it can be 
represented with arbitrary accuracy by a string of bits. 
However, there are situations where information cannot 
be encoded in such a way. For example, the emitter (con- 
ventionally called Alice) wants to indicate to the receiver 
(Bob) a direction in space. If they have a common co- 
ordinate system to which they can refer, or if they can 
create one by observing distant fixed stars, Alice simply 
communicates to Bob the components of a unit vector n 
along that direction, or its spherical coordinates 9 and 
4>. But if no common coordinate system has been estab- 
lished, all she can do is to send a real physical object, 
such as a gyroscope, whose orientation is deemed stable. 

In the quantum world, the role of the gyroscope is 
played by a system with large spin. For example, Al- 
ice can send angular momentum eigenstates satisfying 
n ■ J\ip) — j\ip)- This is essentially the solution pro- 
posed by Massar and Popescu jij who took N parallel 
spins, polarized along n. The fidelity of the transmission 
is usually defined as 

F=(cos 2 (x/2)) = (l + (cos X ))/2, (1) 

where x is the angle between the true n and the direction 
indicated by Bob's measurement. The physical meaning 
of F is that 1 — F = (sin 2 (%/2)) is the mean square error 
of the measurement, if the error is defined as sin(x/2). 
The experimenter's aim, minimizing the mean square er- 
ror, is the same as maximizing fidelity. We can of course 
define "error" in a different way, and then fidelity be- 
comes a different function of x an d optimization leads to 
different results E). Here, we shall take Eq. (1) as the 
definition of fidelity. 

Massar and Popescu showed that for parallel spins, 
1 — F = l/(N + 2). It then came as a surprise that for 
N = 2, parallel spins were not the optimal signal, and a 
slightly higher fidelity resulted from the use of opposite 
spins ||. The intuitive reason given for this result was 
the use of a larger Hilbert space (four dimensions instead 



of three) . This raises the question what is the most ef- 
ficient signal state for N spins, whose Hilbert space has 
2 N dimensions. Will F approach 1 exponentially? In this 
Letter, we show that the optimal result is a quadratic ap- 
proach, as illustrated in Fig. 1. 

Our first task is to devise Bob's measuring 
method, whose mathematical representation is a pos- 
itive operator- valued measure (POVM) Q. For any 
unit vector n, not necessarily Alice's direction, let 
\j,m(n)) = \j,m(0,4>)) denote the coherent angular mo- 
mentum state JH that satisfies 

J 2 \j,m{n))=j(j + l)\j,m(n)), (2) 

and 

n • J\j,m(n)) = m\j,m(n)). (3) 
We then have [|| 

(2j + l) J de<f,\j,m(e,<f > ))(j,m(6,<t>)\ = l j , (4) 

where 

de<j> ■= sin 6 dO d<j)/4:iT, (5) 

and lj is the projection operator over the (2j + 1)- 
dimensional subspace spanned by the vectors \ j, m(0, </>)). 
If N = 2, so that j is or 1, the two resulting subspaces 
span the whole 4-dimensional Hilbert space. For higher 
N, all the rotation group representations with j < N/2 
occur more than once. We then have, if we take each j 
only once, from or I to N/2, 

Ew+1)= <^±2! _ g±iM + »> , ( o) 

for even or odd N, respectively. For large N, the dimen- 
sionality of the accessible Hilbert space tends to N 2 /4, 
and this appears to be the reason that the optimal result 
for 1 — F is quadratic in N , not exponential. An intu- 
itive argument for this quadratic behavior was given by 
Aharonov and Popescu || . No improvement results if we 
endow the particles with internal quantum numbers such 
as charge or strangeness, so that the entire Hilbert space 
can be spanned by states with distinguishable properties, 
because any additional information that Alice could send 
to Bob would refer to these new quantum numbers, not 
to the direction of n. 

We now turn to the construction of Bob's POVM [Q. 
Let p denote the initial state of the physical system that 
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is measured. All these input states span a subspace of 
Hilbert space. Let 1 denote the projection operator on 
that subspace. A POVM is a set of positive operators 
Efj, which sum up to 1. The index fj, is just a label for 
the outcome of the measuring process. The probability 
of outcome [i is tr(p£J^). In the present case, \i stands for 
the pair of angles 9<f> that are indicated by Bob's measure- 
ment. If we want a high accuracy, these output angles 
should have many different values, spread over the unit 
sphere [Q. For example, the components of a continuous 
POVM, as in Eq. (B), are given by 



E H = (2j + 1) d e<i> \j, m(6, </>))(j, m(9, . 



(7) 



Such a POVM with m = j corresponds to the method 
of Ref. [Q. The choice m = j is not optimal. As shown 
in J! for the case N = 2, signal states with opposite spins 
give a higher fidelity. With our present notations, these 
states are (|0, 0) + 1 1, 0(n)))/ y/2. They involve two values 
of j, but a single value of m, namely 0. 

One possibility to include several values of j in a 
POVM is to take a sum of expressions like (g). This 
brings no advantage, because a convex combination of 
POVMs cannot yield more information than the best one 
of them [|| . Optimal POVM components can always be 
assumed to have rank one. Therefore each one of them 
should include all relevant j: 



where 
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1 \j,m(6,4>)). 



(8) 



(9) 



To verify that this is indeed a POVM, we 
note that in fEg<f, there are diagonal terms (2j + 
l)\j,m(9,<j>)){j,m(6,<f>)\, which give 1 3 , owing to Eq. (g). 
The off-diagonal terms with j\ ^ j% vanish, as can be 
seen by taking their matrix elements between (ji , mi | 
and |j2,TO2) m the standard basis where J z is diagonal. 
We have fl 



(10) 



with a similar (complex conjugate) expression for 
(iii rn(0, 4>)). The rotation matrices V are explic- 

itly given by 



^ a (W) = e <TO *$l(0)e <ma *> 



(11) 



where the Euler angle ip is related to an arbitrary phase 
which is implicit in the definition of |j, m(#, </>)). It is 
crucial that a single value of m occurs in all the com- 
ponents of the vectors \0,(f>) in Eq. (^), so that the un- 
defined phases e ±lm ^ mutually cancel. It then follows 
from Eq. (4.6.1) of Ref. M that all the off-diagonal ma- 
trix elements of J Eg^ vanish, so that we indeed have a 
POVM. 



While Bob's optimal POVM is essentially unique in the 
Hilbert space that we have chosen, Alice's signal state, 
which is 



N/2 
j=m 



(12) 



contains unknown coefficients Cj. The latter are normal- 
ized, 

N/2 

EM 2 = 1 ' ( 13 ) 

but still have to be optimized. 

The probability of detection of the pair of angles 9(f), 
indicated by the POVM component Eg^, is 



(A\E g<j> \A) = d 04 , 



N/2 



Cj^/2j + l{j,m(6,(f>)\j,m(n)) 



We have § 

(j,m(9,cl>)\j,m(n)) = ^ d^ m (x), 



(14) 



(15) 



where x is ^ e angle between the directions n and 9<f), 
and the phase e l7? is related to the arbitrary phases which 
are implicit in the definitions of the state vectors in (15). 
The important point is that e lJ] does not depend on j and 
therefore is eliminated when we take the absolute value 
of the sum in Eq. (UJ). Explicitly, we have 



d® m (x) = cos 2m ( X /2) PjW(co SX ), 



(16) 



where P„ (x) is a Jacobi polynomial We shall 

write x = cos x for brevity, so that the fidelity is 



F=(l + (x))/2. 



(17) 



Our problem is to find the coefficients cj that maximize 
(x) . Owing to rotational symmetry, we can assume that 
Alice's direction n points toward the z-axis, so that dg<$, 
can be replaced by dx/2 after having performed the in- 
tegration over <t>. We thus obtain 



xdx 



N/2 
j=m 



Ci\f2j 



1 



P 



(0,2r. 



\x) 



(18) 



This integral can be evaluated explicitly by using the 
orthogonality and recurrence relations for Jacobi polyno- 
mials @[y] • The result is 
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(19) 
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where Ajk is a real symmetric matrix, whose only non- 
vanishing elements are 

A jj = m 2 /\j{j + 1)1 (20) 

and 

= Aj-U = (f - m 2 )/j^ 2 - 1. (21) 

The optimal coefficients Cj are the components of the 
eigenvector of Ajk that corresponds to the largest eigen- 
value, and the latter is (x) itself. The result of the cal- 
culation is displayed in Fig. 1 for m = (which is best) 
and m = j (which is the method investigated in Ref. jD]). 
For to = and large N, we find that 

1 5.78317/(7V + 3) 2 . (22) 

This ought to be compared to the result of Q , which was 
l/(N + 2). For N — 2 and m = 0, our result coincides 
with Ref. [§. For N = 3, we obtain F = 0.84495 with 
c 3/2 = 0.60362 and Ci/ 2 = 0.79755. The results for larger 
N and intermediate values of to gradually fall between 
those displayed in Fig. 1. Had we chosen a definition of 
fidelity other than Eq. (1), these results would of course 
be different, but the method for solving the problem is 
in principle the same. 

It thus appears that it is advantageous to take the 
lowest possible to (namely to — for even N and m = 5 
for odd N). This is intuitively quite plausible ||. It 
would be interesting and instructive to find a direct proof 
of Eq. (22) that does not rely on a numerical analysis as 
in the present work. 
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FIG. 1. (1 — F) as a function of N. Open circles are for 
m = j (Ref. |l|), closed circles are for to = (this work). 
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